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Abstract A central model in theoretical ecology considers the competition of a range 
of species for a broad spectrum of resources. Recent studies have shown that essen- 
tially two different outcomes are possible. Either the species surviving competition are 
more or less uniformly distributed over the resource spectrum, or their distribution is 
'lumped' (or 'clumped'), consisting of clusters of species with similar resource use that 
are separated by gaps in resource space. Which of these outcomes will occur crucially 
depends on the competition kernel, which reflects the shape of the resource utilization 
pattern of the competing species. Most models considered in the literature assume a 
Gaussian competition kernel. This is unfortunate, since predictions based on such a 
Gaussian assumption are not robust. In fact, Gaussian kernels are a border case sce- 
nario, and slight deviations from this function can lead to either uniform or lumped 
species distributions. Here we illustrate the non-robustness of the Gaussian assumption 
by simulating different implementations of the standard competition model with con- 
stant carrying capacity. In this scenario, lumped species distributions can come about 
by secondary ecological or evolutionary mechanisms or by details of the numerical im- 
plementation of the model. We analyze the origin of this sensitivity and discuss it in 
the context of recent applications of the model. 
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Introduction 



A central model behind th e theoretical description of co mpetition among dissimilar 
species was introduced by iMacArthur and Levins! 1 19671 ). In the model, species are 
characterized by their niche position x±, which measures a trait being relevant for the 
exploitation of a distributed resource. As examples, the niche value X{ may represent 
body size of predators, where the distributed resource are preys with their size distri- 
bution, or Xi could be beak sizes of birds, in which case the resource could be seeds of 
different sizes. Mathematically this leads to a Lotka-Volterra type of competition equa- 
tion, where the competition coefficients are function of the distance between species 
on the niche axis x. This competition kernel is usually taken to be a Gaussian function 
of the niche difference (also called normal curve). The implication of this choice is the 
central topic of this paper. 

The model was originally proposed as part of the hypothesis of limiting similarity, 
namely tha t competing species can coex i st only if they are sufficiently different from 
each other I MacArthur and Levinj|l967l ; lAbrams A mathematical analysis of 

the model revealed that arbitrarily similar species could in fact coexist in som e cases . 
Howev er adding furth er effects to the model, li ke noise dMav a ndMa cArthuij ( 19721 ). 



but see Turclli 



g lurcher ertects to the model, like noise (Mav and Mac 
HjjEl)) or extinction thresholds l|Pigolotti et al. Il2007h 7 



impose a limit 



to the similarity between species. This sensitivity to second order effects has led to the 
conclusion that the model, in its original form , is structurally unstable when used to 
predict limits of similarity I Meszena et aL 20061). The competition model has also been 



applied to describe coevolving species l|MacArthur and Levins! 1 19671; ICase Ill98ll) and 



used in some formulations of the theory of island biogeography 1 Roughgarden 19791 ) 



More recently the same type of model has been si mulated numerically and used a s a ba- 
sis for dynamical models of sy mpatric speciation jDoebeli an d Dieckmann ||2000|) , food 
web assembly and ev olution ijLoeuille and Loreau 1 120051 : 1 Johansson and Ripalbood : 
iLewis and Law elucidating the relation between competition and predator-prey 

interaction s dChesson and KuanglbOQSl) . and for explaining lumped size distributions 
of species ijScheffer and van Nesl 120061 ). For a more e xtensive review of the biol ogical 
applications and the generalization of the model see ( Szabo and Meszena I [20061 ). The 
competition model has been fundamental for the development of basic principles in 
theoretical ecology and it is relevant to achieve a full understanding of its technical 
aspects. 

In almost all applications of the model the chosen competition kernel is Gaus- 
sian. This choice facilitates mathematical analysis, and was justified because the exact 
shape of the kernel was thought to have no influence on the fundamental results of the 
model. However, recent work has shown that the equilibrium solution can be one of 
two fundamentally di fferent types, depending on the form of the competition kernel 
( Pigolotti et al. II2OO7I) . One class of competition kernels preserves all species initially 
introduced in the system, with adjustments only in their relative abundance. The final 
equilibrium is a state with species closely spaced and with roughly similar abundances 
(if the carrying capacity is also uniform) . Another class of competition kernels leads to 
the species being lumped in dense groups (in some cases groups are formed by single 
species), separated by empty regions on the niche axis. Subsequent invasion of new 
species in these 'exclusion zones' is not possible due to competitive exclusion. The 
condition for uniform distribution of species is to have a positive definite competition 
kernel (see definition below) . This criterion is automatically fulfilled wh en the kernel is 
constructed from the overlap of the species utilization of the resource (Roughga rden I 
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1979). If the kernel is not positive definite, a lumpy species distribution with exclusion 
zones emerges. The concern about this discovery is that, even though the Gaussian 
kernel is ecologically sound, it is exactly marginal between the two regimes. This indi- 
cates that numerical inaccuracies and/or secondary ecological effects may violate the 
positive definiteness of the competition kernel and cause a transition from a uniform 
to a lumpy species distribution. 

The objective of this paper is to raise awareness in the theoretical ecology com- 
munity of the potential pitfalls and subtleties associated with the use of Gaussian 
competition kernels or other marginal choices. To illustrate this, the consequences of 
the marginal nature of the Gaussian kernel in the competition model are explored, in 
the idealized case of a uniform carrying capacity. First, the sensitivity to ecologically 
relevant effects that may lead to lumpy distributions are examined. Then we exam- 
ine the sensitivity to the details of numerical implementation. In the last section, we 
discuss the relevance of our results for the applications of the model. 



Model 

The competition model considers n interacting populations, each utilizing a common 
distributed resource x according to a utilization function Ui(x), i = 1, n. The dynam- 
ics of the abundance of species i, Ni, is described by a Lotka-Volterra set of competition 
equations: 

Ni = Ni^l - G v N n .«' = !. •••> n ' C 1 ) 

where the growth rate (considered to be the same for all species) is set to one for 
simplicity, and the carrying capacity K is uniform. Competition in ([1]) is described 
by competition coefficients Gg which are constructed from the overlap of utiliza tion 
functions of competing species ( MacArthur and Levinslll967l ; iRoughgarden Ill979l ): 



J Ui(x)uj(x) dx 
J uf(x) dx 



Gij = J r\l 3 _ s \_ ■ (2) 



A justificati on of (|2|) rests upon considering the probability that consumer i meets 
lsj] 



consumer j |LevinsTl96i ; Roughga rden |[T979l ) 



Often, utilization functions are ignored, and the competition coefficients are pos- 
tulated directly. It is usually assumed that species i has an optimal exploitation of the 
resource at a value x = Xi, and the competition coefficients are taken to depend on the 
difference between the optimal resource values of two competing species, y — \x± — Xj\, 
such that we can introduce the so-called competition kernel, Gij = G(y). Here we use 
a family of competition functions described by a parameter p: 

G ij =G(y)=e-^- x ^", (3) 

which contains the Gaussian kernel when p — 2, or the exponential one when p = 1. 
The width of the kernel a gives the range of competition on the niche axis. Incidentally 
the Gaussian kernel is obtained from Eq. (f2Jl when the utilization functions are also 
Gaussian and of the form Mj = exp(— ((x — Xi)/s) 2 ) with s 2 = er 2 /2. When p < 2 the 
kernels are more peaked around y ~ and for p > 2 they become more box-like (see 
Fig. LI). 
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Note that when competition coefficients are constructed by the formula (|2jl. i.e. from 
the overlap of two utilization functions, the y are always positive definite, meaning that 
a iGijCij > for any set of numbers iRoughgarden Ill979l ) or, equivalently, that 



the Fourier transform of G(y), defined as G(k) = J_°° dx G(x) exp(ikx), does not 
take negative values. This property holds for the family of kernels (|3J for p < 2, but 
not for p > 2 (Fig. The Gaussian kernel is therefore marginal in the sense that, 
corresponding to the limit case p — 2, even a very small perturbation may violate its 
positive definite character, generally believed to be an ecological requirement arising 
from expression ([2}. 

An intuitive explanation for the appearance of the exclusion zones for p > 2 is 
the following. Interaction kernels with large p have a box-like shape. In these cases 
species compete very strongly with other species, roughly within a distance ±a from 
their own niche value. Species with a niche x in that range will therefore not be able 
to invade the resident species, leading to the exclusion zones between them. When p is 
decreased, the resident species compete less and less with neighbouring species, until 
the exclusion zones disappear, leading to the possibility of continuous coexistence. 

Understanding the fact that the transition occurs at p — 2, and also the coexistence 
of more than one species in each cluster, requires a mathematical stability analysis of 
the model. Consider the uniform solution, in which many species having the same 
abundance are densely packed in niche space. Now perturb each population by a small 
quantity ANi, which can be either positive or negative. If the competition kernel 
is not positive defined, there are sets of perturbations such that ANjGjjANj is 
less than zero. One ca n show that such perturbations are amplified by the dynamics 
( Pigolotti et al. 1120071 ). making the uniform solution unstable. The system will then 



evolve to a clustered state, where the distance between clusters is proportional to the 
interaction range a. 

We mention here that a possible generalization is to consider multi-dimensional 
niche spaces. This possibility would complicate the mathematical notation but does 
not introduce qualitative changes. Stability in a multi-dimensional niche space still 
depends on the positive definiteness of the competition kernel. In particular, a multi- 
dimensional Gaussian competition kernel is still marginal in the sense described above. 

We simulated the model JTJ with competition kernel Q for 1000 generations and 
200 species initially at random niche positions. The width of the kernel is a — 0.3 and 
the carrying capacity is K = 10. The niche range is taken to be a; £ [0, 1]. The standard 
mathematical way to avoid effect s due to the borders of the niche space is to adopt 
periodic boundary conditions (e.g. IScheffer and van Nesl (|2006h ). These are introduced 



for mathematical convenience and aim at modeling species far from endpoints in a large 
niche space. Adopting periodic boundary conditions means that niche space is circular, 
so that when the interaction kernels extends beyond the left edge at x = 0, it enters 
back into the right side at x = 1 and vice versa. Periodic boundaries therefore mimic 
an infinite system by considering the niche segment [0, 1] as embedded in an array of 
repeated copies of itself. Mathematically, this is properly implemented by making a 
'kernel wrap', i.e. substitute G(y) in ([3| with G p (y) = ^2 n G(y — n), where the sum 
runs from n — 0, ±1, ±2, ... ± oo. We call this implementation "fully periodic boundary 
conditions", to distinguish it from another possibility considered below. Under fully 
periodic boundary conditions, the stability properties of the uniform solution are the 
same as in the infinite system. 
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Fig. 1 Three interaction kernels (top), their Fourier transform (middle row), and species 
distributions arising from simulation of the model after 1000 generations (bottom). Column 
are (Left) Exponential competition (p = 1); (Center) Gaussian competition (p = 2), and 
(Right) quartic competition (p = 4). Simulations are initiated with 200 species randomly 
distributed, K = 10, and a = 0.3. 



Results 

Simulations using the competition kernel ([3]) with p = 1 (exponential), 2 (Gaussian) 
and 4 (box-like) illustrate the uniform species distributions for p = 1 and p — 2, 
and the lumped species clusters for p = 4 (Fig. [1} . The configurations in Fig. [T] are 
still transient states and, at longer times, configurations with p < 2 become more 
uniform, whereas the periodically spaced clusters of species for p > 2 become thinner 
until they contain only a single species. In any case, from the initial stages until the 
final equilibrium, the main difference between the dynamics for the two classes of 
competition kernel is unchanged: for p < 2 all initial species are preserved, leading 
to dense and evenly distributed configurations, whereas 'exclusion zones' develop for 
p > 2 leading to lumped species distributions. 



Effects of secondary ecological processes 

A natural question is whether the marginal nature of Gaussian competition can be 
brought on by secondary ecological effects. We have checked that adding a small im- 
migration rate does not produce lumpy distributions. Adding noise or an extinction 
threshold (i.e. species are removed when their populations fall below a threshold) result 
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Fig. 2 Final populations after 500 000 generations with speciation and extinction. Species 
whose population goes below 0.1 are removed from the system. Every 100 generations new 
species are introduced close to an existing one. The parent species is chosen with a probability 
proportional to its population; the distance of the new species to its parent is drawn from a 
Gaussian distribution of zero mean and spread a p = 0.02. The new species j is introduced 
with a population uniformly drawn from the interval N £ [2,3]. (Left panel) Gaussian kernel 
(p = 2) and (Right panel) exponential kernel (p = 1). Simulations are performed with fully 
periodic boundary conditions. K = 10 and a = 0.3. 



in a l imit to similarity between species, but does not lead to clustering ( Pigolotti et al. I 
120071 ). This also happens in non marginal cases with p < 2, where the minimum dis- 
tance between species is unrelated to the competition range a. 

Effect of species extinction, invasion, and speciation was simulated by eliminating 
species below a given population threshold, and introducing invading species at a fixed 
rate. If they are introduced at random locations in niche space no patterns are observed. 
If invading speci es are introduced close to e xisting ones, modeling evolutionary change 
and speciation l Lawson and Jensen I [2007). the system ends with a lumped species 
distribution, even for p = 2 (Fig. |2j. However, the same mechanism has no effect 
if an exponential competition kernel (p — 1) is chosen. The interpretation is that 
evolutionary effects favor the formation of lumpy species distributions, but only when 
the competition kernel is close to the Gaussian limiting case. 



Effects due to truncation 



The most obvious numerical simplification is to only partially implement the periodic 
boundary conditions, by omitting the kernel wrap around the niche interval, that is, 
using G(y), with y being the minimum of the two possible distances among species i 
and j (\x.i — Xj\ and 1 — \xj, — Xj\), instead of the periodic kernel G p (y). The resulting 
effective kernel is Gaussian but truncated at \y\ = 1/2, making it no longer positive 
definite. Although the shapes of G(y) and G p (y) are still very similar for the parame- 
ters used here (a = 0.3), the change immediately leads to lumped species distributions 
(Fig. O. In contrast, for p = 1 (or any other values of p < 2 which we have checked), 
changing G p (y) by G(y) has no noticeable effect. Qualitatively, the dynamics for trun- 
cated Gaussian kernels resembles the outcome when the exponent of the competition 
kernel is perturbed just slightly. E.g. using p — 2.1 instead of p — 2 also leads to lumped 



7 



N - 

0.18 - 
0.16 - 
0.14 - 
0.12 
0.1 
0.08 | 
0.06 | 
0.04 [ 
0.02 \ 



Fig. 3 Populations of 200 species after 1000 generations with (left panel) Gaussian competi- 
tion kernel with fully periodic boundary conditions, and (right panel) with truncated Gaussian 
competition kernel (see text). K = 10 and a = 0.3. 



species distributions, even when fully periodic boundary conditions are implemented 
(not shown). 

While the mathematical analysis of "evolutionary diffusion" is rather complicated 
and r esults can be obtained only in the framework of simple models ( Lawson and Jensen I 
2007), the effect of truncation on the different kernels may be analyzed in more depth 
and is discussed in the next subsection. 



Wavelength of Gaussian and exponential instabilities 



We have numerically shown in the previous section that truncation leads to radically 
different results for the Gaussian and exponential kernel. This may be surprising when 
realizing that neither the truncated Gaussian nor the truncated exponential kernel are 
positive defined. The explanation of the different result comes from the wavelength of 
the modes k for which the Fourier transform of the kernel G(k) takes negative value s, 
as it determines the distance between lumps i Pieolotti et al. 1120071 : librt et al. 1120091 ). 

To find these wavelengths, we start from the Gaussian case and show what is the 
effect of truncating the kernel at a distance a (which we take now to be the niche-space 
size, previously scaled to be 1) : 



G p=2 (k) = 



2 2 — k 2g2 

exp(— x /a — ikx) dx — — e 4 



. / a ika 

:rf — I 

\a 2 



+ C.C. 



, (4) 



where erf (a;) = (2/^tt) exp~t 2 dt is the error function, and c.c. stands for complex 
conjugate. We impose a/a 2> 1, so we can get a simpler expression by expanding the 
error function, yielding 



/7TO" 



2e" 



a cos(ka) 2 sin(fca) 
Tva ivka 



(5) 



From the previous expression we see that the wavenumbers that give negative modes 
are order k ~ 2a /a 2 . 
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We now check the same effect in the exponential case: 

/ a 2cr (f — e~~ cos(afc) + e~ ~ ka sin(ak)) 
exp( — \x\ la — ikx) dx = -=— g (6) 

In this case the instability occurs for large wavenumbers, proportional to exp(a/a). 
We checked numerically that kernels with p < 2 behave like the p — 1 case, with the 
wavenumber of the first negative mode growing exponentially with the size of the niche 
space. 

Summarizing, the kernels we considered develop an instability due to the trunca- 
tion, but at very different frequencies: very high for the exponential kernel (exponential 
in the ratio between the system size and the kernel range) compared with the Gaussian 
(proportional to the same ratio). The consequence is that, already when the size of the 
niche space is large but not extremely large compared with the competition distance 
(like the cases considered here, a/a = 0.3), the unstable mode in the exponential case 
has a wavelength being much smaller than the interspecies distance, and is thus unable 
to generate patterns. 



Discussion 



The model Q-© provides a very abstract representation of competition. Both em- 
pirical observations and theoretical approaches, based on explicit consideration of the 
coupled consumer-resou rce dynamics, lead to competition coefficients which are quit e 
different from Gaussian (|Schoenedll974l : IWilson1ll975l : lAckermann and Doebeii"ll2004l ) . 
except in a few particular cases. Even so, the qualitative outcome of the model does 
not depend on the exact shape of the competition kernel, but only on its positive or 
non-positive definite character. We have restricted our considerations primarily to the 
basic model flj with Gaussian interaction kernel and constant carrying capacity since it 
is the simplest implementation, allowing to illustrate in a clear setting the importance 
of G and the issues caused by the choice of a marginal competition kernel. 

The basic model with competition coefficients obtained from the overlap of utiliza- 
tion functions, which yields always positive definite kernels, allows for dense species 
distributions with no limits to similarity. This fundamental solution may be changed by 
three different effects: 1) effects stemming from the competition kernel being no longer 
positive definite lead to lumpy species distributions. Clusters of species will appear, sep- 
arated by exclusion zones in niche space with a spacing proportional to the width of the 
competition kernel a; 2) second order ecological effects like noise, species heterogeneity, 
evolutionary effects, or the introduction of an extinction threshold lead to a limit to 
the similarity with the spacing between species being independent of a; 3) under a non 
constant carrying capacity, patterns of unev enly spaced species, lumpy or not, may 
appear. This lead lSzabo and Meszenal ( 20061 ) to conclude that "the not-very-smooth 
nature of the carrying capacity seems to be essential for limiting similarity". Notice 
that also some types of smooth, but non-uniform, c arrying capacities may originate 
lumped distributions (|Hernandez-Garci'a et al.~ll2000h 

The first case arises when the competition kernel is not positive definite. The main 
point of this paper is that, when the competition kernel is the marginal Gaussian, 
the model becomes very sensitive to additional effects that may lead to lumpy species 
distribution. As an example, we demonstrate d that a simple representation of evo- 
lutionary diffusion ( Lawson and Jensen I [20071 ) may lead to patterns in the Gaussian 
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case. This effect is similar to that of evolutionary d ynamics under assort at ive mating 
which shown to lead to lumpy species distributions ( Doebeli et al. ll2007l ). It is worth 
mentioning that, in the context of evolutionary dynamics, the fact that the presence 
of clust ers depends of the c hoice of the competition kernel has been recently demon- 
strated i Leimar et al. I E0O8V Patterns can also result from numerical approximations, 
such as truncating the tails of a Gaussian competition kernel. This effect is proba- 
bly th e underlying mechanism be hind species clustering observed in recent numerical 
work (Schcffcr and van Nesll2006h . which was used to explain observed lumpy distri- 
butions ( May_et_aL 20071). Additional ecological effects have been identified however 
(Schcffcr and van Nel2 006l which make the species groups a robust phenomenon (See 
iHernandez-Garcfa et al. I (|2000l ) for analytical solutions of this type). In any case, these 
spurious effects can be avoided by paying attention to numerical details or by using a 
competition kernel which is not marginal, e.g. one with p — 1.5, which in practice is 
almost indistinguishable from the Gaussian one. It is worth mentioning that analytical 
(i.e. not numerical) results are not affected by the marginal nature of the Gaussian 
kernel, both in relation to limiting similarity ( May and MacArthui 19721). coevo lution 
( Case |[l98ll ) or criteria for sympatric speciation ( Doebeli and Dieckmannll2000l ). The 
marginal nature of Gaussian com petition kernel may however affect numerical work on 
food web evolution and asse mbly ( Doebeli and Dieckmann~ll2000l ; lLoeuille and Loreau I 
l2005l ; lLewis and Law llioOTl ). 

Since a non-positive definite competition kernel leads to lumpy species distributions 
a natural question is whether a non-positive definite kernel can result from hypothesis 
on ecological interactions. This case is often neglected in the literature, since assuming 
Eq. Q automatically leads to a positive definit e com petition kernel (|Roughgarden I 
1979). However, as emphasized in lMeszena et all (|2006l ) and references therein, under 
quite general assumptions one should introduce two different utilization-like functions: 
a sensitivity function Si(x), describing the effect of the resource at x on the growth of 
species i, and an impact function Di(x), describing the depletion of resources produced 
by i. Then, the competition coefficients depend on the overlap of these two quantities 
J Si(x)Dj(x)dx, and reduce to ([2]) only if the sensitivity and impact functions are 
proportional, with the constant of proportionality being the ecological efficiency. When 
the ecological efficiency is a function of x, and the sensitivity and impact functions 
are no longer proportional, the competition kernel may cease to be positive definite 
( Hernandez-Garcfa et al~ll2000l ). 

The third m echanism is that of a non-co nstant carrying capacity K(x), which has 
been explored bv lSzabo and Meszena I (|2006l ). They found that some choices of carrying 
capacity leads to an irregular species lumping. The effect of non-constant carrying 
capacity in conju nction with both positive and non -positive definite competition kernels 
was explored by IHernandez-Garcfa et al.l ( 2000l) . The emerging picture is that the 
two mechanisms are independent. The cases in which a non-constant carrying capacity 
leads to uniform species distributions can also be destabilized by a non-positive definite 
kernel. This means that the mechanism explored here is not a particularity of constant 
carrying capacity but is present also in more general settings. It has been suggested 
that the effect of cutting the ke rnel is similar to that of having a finite niche space, i.e. 
of a box- like carrying capacity dFort et al. 1120091 ). However, this analogy does not hold 
to a rigorous mathematical analysis, which has demonstrated that lumps may appear 
when the carrying capacity is non-constant, but that this e ffect is not mathematically 
equivalent to a cut ker nel ( Hernandez-Garcfa et al.~ll2000h . Moreover, recent results 
( Baptestini et al.ll2009l ) show how the effect of closed boundaries depend on the details 
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of the boundary shape, so that the box-like case is not representative of a general finite 
niche space. 

Having outlined the reasons that may cause the three different outcomes, the ques- 
tion arises if it is possible to infer whether one effect or the other is at play from the 
result of a numerical integration of the competition model. It can be difficult to distin- 
guish between a uniform distribution of discrete species and a lumpy one when lumps 
are very narrow and close. Here, the fact that in the lumpy distribution the spacing 
of the lumps is proportional to the width of the competition kernel a can be used. 
If changing a results in a change in the distance between species proportional to a, 
the effect is due to a non-positive definite competition kernel and vice versa. In the 
case where the effect is due to th e carrying capacity being n on-constant the spacing of 



species is usually more irregular 1 Szabo and Meszena 112009 ) 



To summarize: in line with previous works we have found that the case of continu- 
ous coexistence (no limits to similarity) may be limited by a variety of effects, specially 
for the Gaussian kernel which has a marginal character. We have underlined that there 
are different ways to limit similarity, some leading to lumpy species distributions and 
others not. We hope that this article will increase the awareness in the theoretical 
ecological community of the potential subtleties associated with the use of the Gaus- 
sian competition kernel. Even though this functional form appears to be natural, in 
particular for analytical work, it may not be the most prudent choice for numerical 
exploration of the niche model. 
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